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ABSTRACT 

Context. The earliest phases of clustered star formation and the origin of the stellar initial mass function (IMF) are currently much 
debated. In one school of thought the IMF of embedded clusters is entirely determined by turbulent fragmentation at the prestellar 
stage of star formation, while in a major alternative view it results from dynamical interactions and competitive accretion at the 
protostellar stage. 

Aims. In an effort to discriminate between these two pictures for the origin of the IMF, we investigated the internal and relative 
motions of starless condensations and protostars previously detected by us in the dust continuum at 1.2 mm in the L1688 protocluster 
of the Ophiuchus molecular cloud complex. The starless condensations have a mass spectrum resembling the IMF and are therefore 
likely representative of the initial stages of star formation in the protocluster. 

Methods. We carried out detailed molecular line observations, including some N2H + (l-0) mapping, of the Ophiuchus protocluster 
condensations using the IRAM 30m telescope. 

Results. We measured subsonic or at most transonic levels of internal turbulence within the condensations, implying virial masses 
which generally agree within a factor of ~ 2 with the masses derived from the 1.2 mm dust continuum. This supports the notion that 
most of the L1688 starless condensations are gravitationally bound and prestellar in nature. We detected the classical spectroscopic 
signature of infall motions in CS(2-1), CS(3-2), IFCO(2i2- In), and/or HCO + (3~2) toward six condensations, and obtained tentative 
infall signatures toward 10 other condensations. In addition, we measured a global one-dimensional velocity dispersion of less than 
0.4 km s _I (or twice the sound speed) between condensations. The small relative velocity dispersion implies that, in general, the 
condensations do not have time to interact with one another before evolving into pre-main sequence objects. 

Conclusions. Our observations support the view that the IMF is partly determined by cloud fragmentation at the prestellar stage. 
Competitive accretion is unlikely to be the dominant mechanism at the protostellar stage in the Ophiuchus protocluster, but it may 
possibly govern the growth of starless, self-gravitating condensations initially produced by gravoturbulent fragmentation toward an 
IMF, Salpeter-like mass spectrum. 

Key words, stars: formation - stars: circumstellar matter - ISM: clouds - ISM: structure - ISM : kinematics and dynamics - ISM: 
individual objects (L1688) 

1. Introduction clouds such as the L1688 clump in Ophiuchus (also known as 

the p Ophiuchi main cloud; d ~ \5Q pc), the Serpens central 
While most stars are believed to form in clusters (e.g. clump (d ~ 300 pc), or the NGC 2068/2071 clumps in Orion B 
l Adams & Myers || 200U | Lada & Lada| | 2003|) , our present theoret- {d „ 400 pc) have uncovere d 'complete' samples of prestel- 
ical understanding of the star formation process is essentially lar condensations whose associated mass distributio ns resem- 
hmited to isolated dense cores and protostars (e.g. |Shu et al.| ble the stellar IMF (e.g. Motte, Andre, & Ne ri 1998 -|MA N98; 
[19871 |2004j). Studying the formation and detailed properties |Testi & Sargenlll998t iJohnstone et al.ll2000l iMotte etal.1 12001 : 
of prestellar condensations in cluster-forming clouds is thus of Bontemps et al. 200l[ and references therein). In particular, this 
prime importance if we are to explain the origin of the stellar is me ca se for the population of 57 starless condensations iden- 
initial mass function (IMF). Some progress has recently been tified by | MAN98 | in their 1 .2 mm continuum mosaic of L1688 
made in this area (cf. * ' "' ' ^' " " """"" 



Mott e & AndrefeOOl 



Andre, Ward-Thompson, & B arsony|2000t with me MPIfR bolometer array (MAMBO) on the IRAM 30m 
| Ward-Thompson et al. || 2007| for reviews). te l e scope. These Ophiuchus condensations, which were identi- 



On the observational side, recent ground-based (sub)- ne d using a multi- resolution analysis equivalent to a wa velet 
millimeter continuum surveys of a few nearby cluster-forming decomposition (cf. IStarck et all 119981: IMotte et all 120031) . are 
seen on the same spatial scales as protostellar envelopes (i.e., 



Send offprint requests to: pandre@cea.fr ~ 2300^-500 AU or ~ 15"-30" in L 1688). Their mass spectrum 
* Based on observations carried out with the IRAM 30-meter tele- is consistent with the Sal peterl(ll955l) power-law IMF at the high- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany), mass end and shows a tentative break at ~ 0.4 M (see Fig. [TJ. 
and IGN (Spain). 
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This break is reminiscent of t he flattening observed in the IM F of 
field stars below 0.5 M (e. g. lKroupal200 U IChabrierl2003l) . also 
present in the mass function of LI 688 pre-ma in sequence objects 
dLuhmanet aT1l2000t iBontemps et~alll2001l) . If real, the break 
occurs at a mass comparable to the typical Jeans mass in the 
dense («H2 ~ 10 5 c m" 3 ) DCO + cores of the central O phiuchus 
molecu lar cloud (cf. lLoren. Wootten & Wil king 1990). The re- 
sults of IMAN98] in L1688 were essentially confirmed by inde- 
pendent (sub)-millimeter dus t continuum surveys of the same re- 
gion with SCUBA on JCMT dJohnstone et alj20"00l) and SIMBA 
on SEST (StankietaL||2006j). 
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Fig. 1. Cumulative mass distribution of the 57 starless condensa- 
tions identified by MAN98 in the L1688 protocluster (histogram 
with error bars). For com parison, the sol id curve shows the shape 
of the field star IMF (e.g. lKroupall200ll) . while the dashed curve 
corre sponds to the IMF of multiple systems (e.g. Chabrier 2003, 
2005). The star markers represent the mass function of (pri- 
mary) pre-main sequence ob jects in L1688 as deriv ed from the 
ISOCAM mid-IR survey of IBontemps et all (|2001|) . The dotted 
line shows a iV(> M) cc M~ 06 power-law distribution corre- 
sponding to the typical mass s pectrum found for CO clumps (see 
iBlidl 19931: iKramer et al.lll998l) . Note the flattening of the mass 
distributions below ~ 0.4 M G and the apparent excess of starless 
condensations over stellar systems at the low-mass end. 



Such a close resemblance of their mass spectrum to the IMF, 
in both shape and mass scale, su ggests that the starless con- 
densations identified by MAN98 are about to form stars on a 
one-to-one or perhaps one-to-two basis, with a high local effi- 
ciency, i.e., M+fMpm i 50%. This strongly supports scenarios 
according to which the bulk of the IMF is at leas t partly deter- 
mined by pre-collapse c loud fragmentation (e.g. iLarsonl 1985, 
120051: lElmegreen|[T997l: lPadoan& Nordlund! 120021) . The prob- 
lem of the origin of the IMF may thus partly reduce to a good 
understanding of the processes responsible for the formation 
and evolution of prestellar condensations. Additional processes 
are likely to be required, however, to account for the forma- 
tion of binary/multiple systems and fully explain the low-mass 
(M < 0.3 M Q ) end of the IMF. Indeed, while m ost young stars 
are ob served to be in close multiple systems (e.g.lDuchene et al] 
120041) . the 1.2mm continuum survey of IM AN981 did not have 
enough spatial resolution to probe multiplicity within the L1688 
condensations. Furthermore, multiple systems are believed to 
form after the prestellar stage by subsequent dynamical frag- 
mentation during th e collapse phase, clos e to the time of proto- 
star formation (e.g. Goodwin et al. 2007). One would thus ex- 



pect the masses of the Ophiuchus prestellar condensations to be 
more directly related to the masses of multiple systems than to 
the masses of individual stars. Surprisingly, the shape of the con- 
densation mass spectrum agrees better with the IMF of individ- 
ual field stars (solid curve in Fig. 1) than with the IMF of mul- 
tiple systems (dashed curve in Fig. 1). Clearly, the link between 
the condensation mass spectrum and the IMF is less robust at the 
low-mass end than at the high-mass end. 

On the theoretical side, two main scenarios have been 
proposed for clustered star formation in turbulent molecular 
clouds. In the first scenario, the distribution of stellar masses 
is primarily determined by gravoturbulent cloud fragmenta- 
tion at the prestellar stage. Briefly, self-gravitating conden- 
sations form as turbulen c e-generated density fluctua tions (e.g. 
iKlessen & Burkertl 120001 IPadoan & Nordlu nd 2002), then de- 
couple from their turbulent environment through the dissipa- 
tion o f MHD waves on small scales (e.g. lNakano| [l998; Mvers 
1998), and eventually collapse with little interaction with their 
surroundings. Such protocluster condensations are local minima 
of turbulence (traced by narrow linewidths) and have small rel- 
ative motions with respect to one another and to the surround- 
ing gas. A given star is entirely formed from (a fraction of) the 
gas that was initially present in the corresponding prestellar con- 
densation. Thus, in this scenario, the IMF derives directly from 
the co ndensation mass distribution (CMD) (Pad oan & Nord lund 
|2002|) . which provides a simple explanation for the obse rved 
similarity between the CMD and the IMF (e.g. IMAN98I and 
Fig.ffll. 

By contrast, in the second scenario, the distribution of stellar 
masses resu lts entirely from the dynam ics of the parent proto- 
cluster (e.g. lBonnell et alj|1998ll2001bl) . Here, a protocluster is 
initially made up of gas and protostellar seeds. These protostel- 
lar seeds result from turbulence-generated cloud structure like 
in the first scenario, but their initial masses are unrelated to final 
stellar masses. The seeds travel in the gravitational potential well 
of the system and are characterized by a large, essentially virial 
velocity dispersion. They accrete mass competitively as they ex- 
ecute several orbits through the protocluster. The seed trajecto- 
ries within the protocluster are highly stochastic in nature and 
feature close encounters, merging and/or dynamical ejections. 
In this alternative scenario, competitive accretion and dynami- 
cal interactions between individual protocluster members play 
a dominant role in shaping the resulting IM F at the protostellar 
(Class 0/Class I) stage (e.g. lBate et alj|2003l) . Furthermore, most 
of the final mass of a given star comes from gas that was ini- 
tially not g ravitationally bound to the co rresponding protostellar 
seed(s) (cf. lBonnell. Vine. & Batell2004l) . 

In an effort to discriminate between these two broad pictures 
for the origin of the IMF and furth er constra in the nature of the 
starless condensations identified by MAN98, we carried out de- 
tailed molecular line observations of the central Ophiuchus pro- 
tocluster with the IRAM 30m telescope. The present paper de- 
scribes the results of these line observations and discusses them 
in the context of the above-mentioned theoretical scenarios for 
clustered star formation. 

The layout of the paper is as follows. Section|2]provides obser- 
vational details. Section [3] presents the results of our line map- 
ping observations. Section|4]analyzes the constraints set by these 
observations on the kinematics of the L1688 protocluster. We 
discuss the implications of our results for our understanding of 
cluster-forming clouds in Sect. [5] Our conclusions are summa- 
rized in Sect. [6] 
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Table 1. Adopted line rest frequencies and telescope efficiencies. 
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Transition 


Frequency" 




Ref. c 


HPBW rf 


¥ € 

eff 


< 


B e ff 




(MHz) 


(km s" 1 ) 




(") 




(2000) 


(1998) 


H 1J CO + (1-0) 


86754.294(30) 


0.10 


(1) 


28.4 


0.92 


0.77 




N 2 H + (101-012) 


93176.265(7) 


0.023 


(2) 


26.4 


0.92 


0.77 


0.73 


C 34 S(2-1) 


96412.952(1) 


0.003 


(3) 


25.5 


0.92 


0.77 


0.73 


CS(2-1) 


97980.953(1) 


0.003 


(3) 


25.1 


0.92 


0.80 


0.73 


H 2 CO(2 12 -l H ) 


140839.518(7) 


0.015 


(1) 


17.5 


0.90 


0.65 


0.54 


DCO + (2-l) 


144077.319(50) 


0.10 


(4) 


17.1 


0.90 


0.65 




CS(3-2) 


146969.026(1) 


0.002 


(3) 


16.7 


0.90 


0.65 




HCO + (3-2) 


267557.625(17) 


0.019 


(1) 


9.2 


0.85 


0.49 





(a) The frequency uncertainty in units of the last digit is given in parentheses. 

(b) Frequency uncertainty converte d to velocity u nits. 

(c) Reference for the frequency: (l) lLovasl dl992l) . (2) ICaselli et"ai] j!995l) . (3) lGottlieb e"t"aH | |2003|) . (4) |Pickett et"ai] jl998h . 

(d) Half -power beamwidth (HPBW) of the IRAM 30m telescope. 

(e) The forward efficiency was 0.95 for N 2 H + (l-0) in 2005 (OTF map in Oph Bl). 
(/) The main-beam efficiency was 0.77 for N 2 H + (l-0) in 2005. 



2. Observations 

We used the IRAM 30m telescope at Pico Veleta, Spain, in June 
1998, July 2000, and June 2005 to carry out millimeter line ob- 
servations of the DCO + cores and 1.2 mm continuum condensa- 
tions of L1688 in the following molecular transitions: 
N 2 H + (l-0), H 13 CO + (1-0), CS(2-1), C 34 S(2-1) at 3 mm, CS(3-2), 
H 2 CO(2 12 -l n), DCO + (2-l) at 2 mm, and HCO + (3-2) at 1.1 mm. 
Our adopted set of rest line frequencies is given in Table Q] The 
half-power beamwidth of the telescope was ~ 26 ", ~ 17 ", 
and ~ 9 " at 3 mm, 2 mm, and 1.1 mm, respectively. We used 
four SIS heterodyne receivers simultaneously and an autocor- 
relation spectrometer as backend, with a spectral resolution of 
20^40 kHz at 3 mm and 2 mm, and 40 kHz at 1 . 1 mm. The cor- 
responding velocity resolution ranged from 0.04 to 0.07 km s 
per channel, depending on the observed transition. All our obser- 
vations were performed in single sideband mode, with sideband 
rejections of 0.01, 0.1 and 0.05 at 3 mm, 2 mm and 1.1 mm, re- 
spectively. The resulting calibration uncertainty is ~ 10%. The 
forward and beam efficiencies of the telescope used to convert 
antenna temperatures T*^ into main beam temperatures r m b are 
listed in Table Q] The telescope pointing was checked every ~ 2 
hours on NRAO 530 and/or 1514-241 and found to be accurate 
to ~ 4" (rms). The telescope focus was optimized on 3C273 and 
NRAO 530 every ~ 2 hours on average. Single-point line obser- 
vations were performed in the position switching mode, while 
extensive mapping was performed in the "on-the-fly" (OTF) 
mode. All of the data were reduced with the CLASS software 
packagf]. 

3. Molecular line mapping results and analysis 

3.1. Detections of protocluster condensations in N 2 H + (1-0) 
and other tracers 

The N2H + (l-0) integrated intensity mapfl we took toward the 
six main DCO + cores of the Ophiuchus central cloud are shown 
in Fig. 12 The positi ons of the starless 1.2 mm continuum con- 
densations found by MAN98 are marked by crosses. N 2 H + (1- 
0) emission was found to be present toward most of these con- 
densations. Altogether, we performed N2H + (l-0) observations 

1 see |http://www.iram.fr/ IRAMFR/ GnTjASl 

2 Our N 2 H + (l-0) data cubes are available in FITS format at the CDS. 



to ward 48 o f the 57 compact starless condensation^ identified 
by IMAN98I and detected line emission for 41 of them, either 
through position-switch integrations (rms ~ 0. 1 K in units) or 
in OTF maps (rms ~ 0.1-0.2 K, depending on core). However, 
the ~ 26" angular resolution of our N 2 H + observations is not 
always sufficient to clearly distinguish between line emission 
from the condensations themselves and emission from the par- 
ent cloud/DCO + core. In order to estimate the background N 2 H + 
emission from the parent core, we used our OTF maps to produce 
a smoothed N 2 H + (l-0) image of the local background toward 
each object. In pra ctice, we used the m ulti-resolution wavelet 
analysis software of IStarck et al.l (Il998l) to decompose our OTF 
N 2 H + (l-0) data cubes into two wavelet "views", corresponding 
to small (~ 20-40") and large ( i 40") spatial scales, respec- 
tively. The background emission was estimated from the "large- 
scale" view, while an estimate of the background-subtracted 
emission from the condensations was provided by the "small- 
scale" view. Such a multi- resolutio n decomposition is similar to 
the analysis performed by |MAN98l to separate compact conden- 
sations and protostellar envelopes seen on angular scales £ 30- 
60" from more extended ( ~ 1') cloud structure in their 1.2 mm 
dust continuum image. The fact that most of the condensations 
detected in OTF maps remained positively detected after back- 
ground subtraction demonstrates that our N 2 H + observations 
successfully probed the condensations themselves and not only 
the dense environment of the parent clump/DCO + cores. 

Additional H 13 CO + (1-0) and DCO + (2-l) OTF maps were 
obtained toward some of the cores (cf. Fig. [3), which were 
background-subtracted and analyzed in a similar fashion. 

At the position of each target condensation/protostar, the 
seven hyperfine components of the N 2 H + (l-0) multiplet were 
fitted simultaneously using the Gaussian HFS (HyperFine 
Structure) fitting routine of the CLASS software package. This 
routine derives the line optical depth by assuming the same ex- 
citation temperature for all hyperfine components, and therefore 
yields an estimate of the intrinsic linewidth (i.e., properly cor- 
rected for optical depth effects provided that the assumption is 
correct). The results of these HFS fits, both before and after 
background subtraction, are given in Table [2] for all the objects 



3 MAN98 counted 58 compact starless condensations but one of 
these - E-MM3 - was subsequently sho wn to be an edge-on T Tauri 
disk as opposed to a prestellar object (e.g. Brandner et al. 2 0001) . 
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Fig. 2. Maps of the N2H + (l-0) intensity, integrated over the seven components of the hyperfine multiplet, taken in the OTF mode 
with the IRAM 30 m telescope toward the DCO + cores Oph B2 (a), Oph Bl (b), Oph A (c), Oph F (d), Oph E (e), and Oph C- 
S (f). The (0,0) offsets correspond to the J2000 equatorial positions (16 h 27 m 27.96,-24°27'06.9"), (16 h 27 m 12 s .41,-24°29'58.0"), 
(16 h 26 m 26 s .45,-24°24'30.8"), (16 h 27 m 24 s .25,-24°40'35.2"), (16 h 27 m 04 s .70,-24°39'12.5"), and (16 h 27 m 01 s .91,-24°34'40.7"), re- 
spectively. The contours go from 2 to 10, 1 to 7, 2 to 16, 1.5 to 6, 0.75 to 3, and 1 to 6 K km s _1 by steps of 2, 1, 2, 1.5, 0.75, 
and 1 K km s _1 , respectively (in T* scale). The angular resolution (HPBW) is shown as a black filled circle. The crosses mark the 
positions of the starless cond ensations identified bvlMAN98l in the dust continuum at 1 .2 mm, the plus symbol the position of the 
N2I-F peak N6 discussed by iDi Francesco et all d2004l) . and the star symbols the positions of Class or Class I protostars. The 
ellipses show the locations of the clumps identified with Gaussclumps in the corresponding N2H + (101-012) background-subtracted 
data cubes (see Sect. 13. H and TableQ. 
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Fig. 3. H 13 CO + (1-0) or DCO + (2-l) integrated intensity maps taken in the OTF mode with the IRAM 30 m telescope toward the 
DCO + cores Oph B2 (a), Oph F (b), and Oph E (c and d). The contours go from 0.3 to 0.9, 0.5 to 1.5, 0.25 to 1, and 0.6 to 2.4 
K km s _1 by steps of 0.3, 0.5, 0.25, 0.6 K km s _1 , respectively (in T* scale). In (a), the negative contours are -0.6 and -0.3 Kkm s . 
To emphasize small-scale structure, large-scale (> 40") emission was subtracted from the H 13 CO + (1-0) map of Oph B2 shown in 
(a) (see Sect. 13. II ). The angul ar resolu tion (HPBW) is shown as a black filled circle. The crosses mark the positions of the starless 
condensations identified by IMAN98I in the dust continuum at 1.2 mm, and the star symbols the positions of Class or Class I 
protostars. The ellipses show the locations of the clumps identified with Gaussclumps in the corresponding background-subtracted 
data cubes (see Sect. 13. II ). 



detected in N2H + (l-0). (Note that, in some cases, the signal-to- 
noise ratio was not good enough to perform a significant HFS 
fit after background subtraction.) The quoted error bars corre- 
spond to standard deviations (lcr) as estimated by the HFS rou- 
tine of CLASS. Examples of N 2 H + (l-0) spectra and HFS fits 
are shown in Fig. [4] A single-component HFS fit failed and a 
two-component HFS fit was required for three condensations 
in Oph Bl, four condensations in Oph B2, one condensation in 
Oph F, as well as the protostar CRBR 85. 

In addition, we a l so use d the Gau ssclumps fitting pro cedure 
of lStutzki & Gustenl dl990t) (see also lKramer et al.ll 19981) to de- 
compose each of our background-subtracted OTF N 2 H + (101- 
012), H 13 CO + (1-0), and DCO + (2-l) data cubes into a series of 
Gaussian-shaped "clumps". The significant (> 5 cr) clumps iden- 
tified with Gaussclumps in the various DCO + cores of L1688 are 
shown as ellipses in Fig.|2]and Fig. [3] Their main characteristics 
are given in Table [3] In this way, a total of 17 starless 1.2 mm 
condensations and 3 Class 0/Class I protostars (VLA 1623, 
LFAM 26, GY 210) were found to have well-defined N 2 H + , 



H 13 CO + , or DCO + counterparts in position-velocity (/— v) space. 
The relevance of such a Gaussclumps decomposition is that it 
further helped us discriminate between the line emission arising 
from the compact condensations themselves (with well-defined 
"positions" in I — v space) and the line emission from the sur- 
rounding, more extended dense gas (with less-well-defined "po- 
sitions" in / - v space). In particular, for six of the above- 
mentioned eight condensations with double N 2 H + HFS com- 
ponents, the Gaussclumps decomposition allowed us to reject 
one of the two components and to identify the component most 
likely associated with the compact 1.2 mm continuum object. In 
Oph Bl, for instance, two HFS components were detected to- 
ward B1-MM2, B1-MM3, and B1-MM4 prior to background 
subtraction (cf. Table [2]), but only one of these two compo- 
nents was found to be associated with each condensation after 
running Gaussclumps on the background-subtracted N 2 H + ( 101 - 
012) data cubes (cf. Table [3]). In the following, we will con- 
sider the multiple components detected toward the remaining 
two condensations with double HFS components (B2-MM2 and 
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Table 2. Results of Gaussian hyperfine fits to the N2H + (l-0) multiplet. 



Prior to background subtraction After background subtraction 



Source 


Coordinates 


cr° 


S/N b 


1 




FWHM 






S/N 6 






FWHM 






^2000 

16 h 


^2000 


(mK) 




(K) 


(kms" 1 ) 


(km s _I ) 




(mK) 




(K) 


(km s~') 


(kms- 1 ) 




(i) 


(2) 


-24° 
(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


A3-MM1 


26:09.7 


23:06 


79 


7 


1.6(2) 


3.243(26) 


0.877(59) 


0.2(24) 














A-MM4 


26:24.1 


21:52 


65 


45 


15.8(1) 


3.243(2) 


0.445(2) 


2.6(1) 


67 


17 


4.5(5) 


3.324(8) 


0.378(18) 


1.5(13) 


A-MM5 


26:25.9 


22:27 


62 


53 


15.3(5) 


3.171(2) 


0.381(6) 


1.9(3) 


64 


15 


9.8(16) 


3.146(5) 


0.183(10) 


10.6(29) 


SM1N 


26:27.3 


23:28 


60 


98 


52.2(8) 


3.573(1) 


0.487(3) 


8.6(2) 


57 


47 


29.0(9) 


3.587(3) 


0.436(6) 


11.4(5) 


SMI 


26:27.5 


23:56 


43 


125 


35.9(0) 


3.645(1) 


0.597(1) 


5.7(0) 


43 


56 


15.2(4) 


3.682(3) 


0.640(7) 


6.1(3) 


A-MM6 


26:27.9 


22:53 


64 


20 


5.0(3) 


3.353(6) 


0.694(16) 


1.4(5) 














SM2 


26:29.5 


24:27 


52 


68 


19.8(3) 


3.509(1) 


0.485(3) 


3.6(2) 


51 


19 


5.8(5) 


3.519(6) 


0.386(13) 


4.9(12) 


A-MM8 


26:33.4 


25:01 


65 


51 


17.3(4) 


3.505(1) 


0.384(4) 


2.7(3) 


64 


16 


3.9(1) 


3.571(6) 


0.332(12) 


0.1(6) 


A-S 


26:43.1 


25:42 


81 


7 


2.9(6) 


3.710(7) 


0.236(20) 


2.2(22) 














VLA1623 


26:26.5 


24:31 


37 


75 


17.1(2) 


3.685(1) 


0.561(3) 


4.4(1) 


37 


30 


6.9(3) 


3.749(4) 


0.532(8) 


4.2(5) 


B1-MM1 


27:08.7 


27:50 


91 


10 


3.1(1) 


4.069(8) 


0.383(21) 


0.1(2) 














B1-MM2 


27:11.8 


29:19 


101 


26 


5.4(5) 


3.387(19) 


0.826(39) 


2.2(7) 


too 


1 1 


1.9(2) 


3.430(42) 


0.981(88) 


0.1(203) 












10.2(8) 


4.050(4) 


0.287(11) 


3.2(9) 






4.8(10) 


4.060(8) 


0.217(26) 


5.1(28) 


B1-MM3 


27:12.4 


29:58 


67 


38 


9.7(8) 


3.287(4) 


0.280(11) 


11.1(15) 
























14.6(4) 


3.821(2) 


0.340(5) 


3.3(4) 


70 


19 


16.9(22) 


3.780(3) 


0.185(7) 


13.7(27) 


B1-MM4 


27:15.7 


30:42 


216 


15 


34.8(68) 


3.551(5) 


0.158(11) 


23.5(54) 


221 


8 


14.8(14) 


3.762(27) 


0.451(37) 


29.8(37) 












14.6(12) 


3.956(5) 


0.341(14) 


2.2(10) 






4.3(7) 


3.967(11) 


0.205(22) 


0.3(13) 


B1B2-MM1 


27:11.3 


27:39 


130 


16 


7.6(7) 


4.070(6) 


0.384(16) 


0.2(12) 


128 


5 


3.0(11) 


4.065(12) 


0.217(39) 


1.5(37) 


BlB2-MM2 f 


27:18.0 


28:48 


126 


6 


2.7(6) 


3.957(25) 


0.556(63) 


2.4(25) 














B2-MM1 


27:17.0 


27:32 


147 


16 


12.6(11) 


4.049(6) 


0.380(14) 


4.5(11) 


144 


4 


6.2(37) 


3.984(15) 


0.184(42) 


15.1(114) 


B2-MM2 


27:20.3 


27:08 


99 


18 


7.5(2) 


3.947(19) 


0.757(34) 


3.7(4) 
























6.7(7) 


4.344(7) 


0.238(18) 


10.3(0) 














B2-MM3 


27:23.7 


28:05 


130 


14 


8.4(2) 


3.785(9) 


0.531(11) 


3.0(1) 
























4.2(9) 


4.376(4) 


0.305(30) 


2.1(23) 














B2-MM4 


27:24.3 


27:45 


262 


14 


21.7(19) 


3.718(7) 


0.393(17) 


6.4(13) 


259 


5 


14.9(46) 


3.678(12) 


0.262(25) 


16.3(66) 












4.6(7) 


4.358(15) 


0.345(38) 


0.1(22) 














B2-MM5 


27:24.9 


27:26 


246 


13 


15.8(4) 


3.693(1) 


0.497(19) 


2.9(0) 


238 


5 


7.9(21) 


3.644(15) 


0.305(28) 


7.7(42) 












4.7(15) 


4.388(62) 


0.597(97) 


8.0(45) 














B2-MM6 


27:25.3 


27:00 


240 


12 


15.6(13) 


3.763(10) 


0.699(26) 


4.0(8) 


248 


4 


10.1(28) 


3.680(22) 


0.424(39) 


14.3(56) 


B2-MM7 


27:27.9 


27:39 


261 


10 


9.2(13) 


4.214(21) 


0.789(60) 


2.0(11) 














B2-MM8 


27:28.0 


27:07 


206 


19 


19.2(11) 


4.144(6) 


0.584(14) 


3.0(6) 


208 


8 


6.3(3) 


4.186(10) 


0.419(17) 


0.1(5) 


B2-MM9 


27:28.8 


26:37 


240 


13 


11.3(12) 


4.093(12) 


0.650(33) 


1.9(10) 














B2-MM10 


27:29.6 


27:42 


162 


15 


8.9(8) 


4.334(9) 


0.569(22) 


0.5(8) 














B2-MM11 


27:29.8 


25:49 


145 


12 


4.2(2) 


4.040(16) 


1.018(51) 


0.1(4) 


146 


4 


4.4(26) 


4.070(54) 


0.633(134) 


15.4(104) 


B2-MM12 


27:29.9 


26:33 


244 


14 


16.4(12) 


4.093(8) 


0.549(20) 


3.5(9) 


255 


4 


8.5(24) 


4.097(12) 


0.238(26) 


7.7(44) 


B2-MM13 


27:32.7 


26:07 


263 


12 


14.1(14) 


3.967(12) 


0.637(29) 


4.1(11) 














B2-MM14 


27:32.8 


26:29 


258 


11 


14.9(14) 


4.162(12) 


0.771(32) 


3.8(9) 














B2-MM15 


27:32.8 


27:03 


267 


15 


31.6(24) 


4.413(5) 


0.356(11) 


7.1(11) 


272 


6 


14.9(36) 


4.445(8) 


0.221(17) 


9.3(41) 


B2-MM16 


27:34.5 


26:12 


55 


54 


15.7(2) 


4.041(2) 


0.626(4) 


3.5(2) 


54 


23 


9.2(7) 


4.074(7) 


0.426(15) 


6.8(10) 


B2-MM17 


27:35.2 


26:21 


250 


12 


13.1(14) 


4.079(12) 


0.622(28) 


3.4(11) 


257 


5 


4.2(14) 


4.110(19) 


0.320(45) 


1.4(33) 


C-W 


26:50.0 


32:49 


154 


12 


14.6(17) 


3.578(4) 


0.192(9) 


7.2(18) 














C-N" 


26:57.2 


31:39 


93 


33 


36.3(16) 


3.811(1) 


0.211(3) 


12.5(8) 














C-MM2 


26:58.4 


33:53 


263 


8 


18.2(27) 


3.879(12) 


0.429(24) 


11.5(26) 














C-MM3 


26:58.9 


34:22 


281 


9 


35.2(46) 


3.910(7) 


0.307(15) 


17.4(30) 














C-MM4 


26:59.4 


34:02 


287 


9 


28.6(49) 


3.852(9) 


0.334(18) 


15.7(38) 














C-MM5 


27:00.1 


34:27 


62 


38 


38.5(12) 


3.809(1) 


0.314(3) 


17.6(7) 


62 


16 


22.1(7) 


3.803(5) 


0.238(6) 


30.0(25) 


C-MM6 


27:01.6 


34:37 


61 


38 


37.1(22) 


3.710(1) 


0.327(6) 


17.4(13) 


61 


18 


24.8(26) 


3.667(4) 


0.249(7) 


26.3(33) 


C-MM7 


27:03.3 


34:22 


294 


6 


34.1(15) 


3.691(10) 


0.294(11) 


30.0(44) 














E-MMI' 


26:57.7 


36:56 


153 


5 


2.1(2) 


4.434(21) 


0.360(39) 


0.1(195) 














E-MM2d 


27:04.9 


39:15 


61 


39 


13.5(4) 


4.496(1) 


0.287(4) 


3.6(4) 


62 


21 


9.0(7) 


4.482(3) 


0.253(8) 


5.5(10) 


E-MM4 


27:10.6 


39:30 


41 


26 


5.1(3) 


4.232(3) 


0.331(7) 


1.7(6) 


41 


16 


4.8(6) 


4.219(6) 


0.254(13) 


6.7(19) 


F-MM1 


27:22.1 


40:02 


59 


26 


8.9(5) 


4.726(5) 


0.476(11) 


4.6(7) 


62 


15 


3.4(1) 


4.789(6) 


0.329(13) 


0.1(9) 


F-MM2 


27:24.3 


40:35 


43 


53 


18.2(2) 


4.133(0) 


0.229(4) 


10.8(0) 


44 


22 


6.4(8) 


4.148(4) 


0.188(11) 


9.6(23) 












7.7(2) 


4.592(2) 


0.354(6) 


0.6(3) 






4.6(4) 


4.584(1) 


0.304(11) 


2.0(10) 


CRBR85 


27:24.7 


41:03 


215 


13 


25.4(26) 


4.067(4) 


0.230(10) 


10.2(17) 


213 


6 


12.3(28) 


4.032(8) 


0.200(19) 


11.1(41) 












3.3(3) 


4.649(20) 


0.402(44) 


0.1(4) 














YLW15 


27:27.0 


40:50 


37 


94 


29.5(4) 


4.162(1) 


0.258(2) 


8.1(2) 


36 


54 


19.0(10) 


4.172(2) 


0.212(4) 


10.1(9) 



Notes: the numbers in parentheses indicate the uncertainty in units of the last digit. 

(a) cr is the rms noise in T* scale. 

(b) S/N is the signal-to-noise ratio. 

(c) The fitting function is T*(v) = ^(1 - e - T{r) ). In the optically thick case, Pi = T*' pmk x T tot . 

(d) Ttot is the total optical depth o f the N2H + (l-0) multiplet. The optical depth of the isolated component 101-012 is 

(e) Composite starless clump (cf. MAN98) not included in the condensation mass distributions shown in Fig.[TJand Fig. [8] 



B2-MM4) as representative of independent objects when we dis- 3.2. Linewidths and virial masses of the protocluster 
cuss the statistics of linewidths (Sect. 13.21 ) and relative motions condensations 
(Sect. 14.31 ) between condensations. 

Based on the measured N2H + (l-0) linewidths (Tabled, we can 
calculate the nonth ermal component of the line-of-sight veloc- 
ity dispersion (cf. Mvers 1999) and estimate a virial mass for 
each condensation. Table [4] gives the deconvolved (FWHM) di- 



Andre et al.: Kinematics of the Ophiuchus protocluster condensations 
Table 3. Compact objects identified with Gaussclumps in the background-subtracted N2H + ( 101 -012) data cubes. 
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Object 


Aa 


AS 


Vo 


a 


Ax 


Ay 


4> 


FWHM 


Av 


Identification 


name 


C) 


(") 


(kms" 1 ) 


(K) 


(") 


(") 


C) 


(") 


(km s" 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


A-l 


6.0 


63.6 


3.52 


2.16 


27.00 


72.50 


81.2 


35.8 


0.49 


A-SM1N 


A-2 


67.1 


-16.6 


3.42 


1.85 


28.40 


44.70 


46.6 


33.9 


0.37 


A-N6 a 


A-3 


-2.6 


144.1 


3.14 


1.01 


27.00 


60.90 


61.5 


34.9 


0.34 


~ A-MM5 


A-6 


10.9 


39.7 


3.90 


0.65 


27.00 


38.70 


112.1 


31.3 


0.22 


A-SM1 


A- 15 


0.4 


8.9 


3.89 


0.35 


27.00 


27.30 


25.1 


27.1 


0.13 


VLA 1623 


Bl-1 


38.9 


-33.0 


4.00 


0.73 


31.80 


32.20 


50.6 


32.0 


0.37 


B1-MM4 


Bl-4 


-0.8 


0.0 


3.75 


0.60 


27.00 


27.40 


86.3 


27.2 


0.25 


B1-MM3 


Bl-20 


-22.3 


53.8 


4.24 


0.45 


27.00 


58.10 


41.8 


34.6 


0.15 


B1-MM2 


B2-1 


67.8 


9.3 


4.43 


1.09 


30.80 


48.20 


104.2 


36.7 


0.29 


B2-MM15 


B2-2 


56.5 


50.5 


3.91 


0.92 


46.40 


27.30 


-6.3 


33.2 


0.35 


B2-MM13 


B2-3 


16.6 


1.0 


4.17 


0.78 


34.40 


50.20 


109.2 


40.1 


0.31 


~ B2-MM8 


B2-4 


72.3 


-39.3 


4.02 


0.77 


30.20 


43.20 


134.7 


35.0 


0.36 




B2-5 


-30.0 


-10.2 


3.76 


0.66 


51.70 


28.60 


-36.5 


35.4 


0.37 


B2-MM5 


B2-6 


49.4 


-33.3 


3.53 


0.54 


27.10 


35.80 


121.9 


30.6 


0.46 




B2-7 


-20.9 


19.9 


3.41 


0.49 


27.00 


46.40 


94.9 


33.0 


0.24 




B2-8 


83.7 


45.4 


4.16 


0.51 


27.00 


33.90 


17.7 


29.9 


0.28 


~ B2-MM16,17 


B2-10 


-70.5 


-49.9 


3.79 


0.45 


27.00 


48.40 


171.4 


33.4 


0.22 


B2-MM3 


B2-11 


10.1 


-36.9 


4.41 


0.45 


27.00 


27.30 


78.7 


27.1 


0.23 


~ B2-MM7.10 


C-l 


-2.4 


-6.2 


3.67 


0.93 


63.70 


32.10 


15.1 


40.5 


0.32 


C-MM6 


C-2 


-39.1 


21.4 


3.92 


0.76 


50.50 


39.80 


11.7 


44.2 


0.31 


C-MM3, C-MM5 


C-3 


-21.5 


72.0 


3.78 


0.42 


27.00 


60.70 


140.6 


34.9 


0.19 




E-l 


7.7 


3.3 


4.52 


0.61 


50.30 


34.00 


5.5 


39.8 


0.25 


E-MM2d 


E-4 


80.6 


-10.9 


4.28 


0.37 


27.00 


27.30 


24.7 


27.2 


0.28 


E-MM4 


F-l 


33.7 


-22.5 


4.15 


0.82 


66.10 


30.60 


6.0 


39.3 


0.24 


~ YLW 15 


F-2 


-8.1 


29.9 


4.54 


0.51 


27.00 


42.10 


60.3 


32.1 


0.26 


~ F-MM1 


F-4 


5.2 


-12.9 


4.66 


0.41 


27.00 


35.30 


62.0 


30.3 


0.31 


~ F-MM2 



Notes: In col [1], the first part of the object name indicates the DCO + core in which the object is embedded. Columns [2] and [3] are the offsets 
of the object with respect to the core position given in the caption of Fig. [2] Columns [4] and [5] are the center velocity and peak intensity, 
respectively. Columns [6] and [7] give the major and minor angular diameters of the two-dimensional Gaussian fitted to the object, and col. [8] 
is the position angle of the minor axis. Columns [91 an d [10] are the typical FWHM size and FWHM line width, respectively, 
(a) N 2 H + peak N6 discussed bv lDi Francesco et al.l i2004h . 



(a) 




^vW^^wV 



B2-MM8 (°) - 

after background subtraction 

5 10 15 20 

Velocity (km s ) 



E-MM2d < d ' 
after background subtraction 



5 10 15 

Velocity (km s 1 



Fig. 4. Examples of N2H + (l-0) spectra and Gaussian HFS fits 
obtained toward the condensations B2-MM8 (left) and E-MM2d 
(right). The spectra shown in the top row are those observed to- 
ward B2-MM8 (a) and E-MM2d (b) prior to background sub- 
traction, while the spectra shown in the bottom row are those 
obtained toward the same sources after subtracting the local 
background emission estimated using a multi -resolution wavelet 
decomposition (see text). 



ameter (col. [2], [MAN98) and gas + du st mass derived from 
the 1.2 mm continuum (col. [41. IMAN98I) . the assumed gas/dust 



temperature (col. [3], cf. MAN98), the thermal velocity disper- 
sion o-jip) for a particle of mean molecular weight p = 2.33 
(col. [5]), the one-dimensional nonthermal velocity dispersion 
cr NT obtained from the N2H + (l-0) linewidth (col. [6]), the non- 
thermal to thermal velocity dispersion ratio cr NT /ctt (col. [7]), 
the estimated virial mass M vi > (col. [8]), and the virial mass ratio 
a V i r = M V j r IM\2mm (col. [9]) for each observed object. Columns 
[10] to [13] list the values of the same parameters when esti- 
mated from the background-subtracted spectra. For some con- 
densations, the dust temperatures (col. [3]) adopted to estimate 
M\2mm (col. [4]) from t he measu red 1.2mm flux differ slightly 
from those assumed by IM AN981 Here, we used Tj - 12 K for 
all condensations except SM1N, SMI, A-MM6, SM2, A-MM8 
in Oph A, which are likely w armer (Tq = 20 K ) due to their 
proximity to the B3 star SI (cf. lAndre et al.llT993b . Our adopted 
default Tj - 12 K is consistent with the typical mass-averaged 
temperatu re found by dust radiative transfer models of star less 
cores (cf. lAndre et"aUl2003t IStamatellos & Whitworfhll2003l) . In 
agreement with the radial profiles measured by |MAN98l the val- 
ues quoted for M v { r and ar vi > assume that the density structure of 
the sources approaches that of centrally-condensed spheres with 

outer profiles such as p cc r 
Here, cr, , 



2 , for which M vir * 3r OTlrf ^. 

2 NT is the total (thermal + nonther- 
mal) line-of-sight velocity dispersion and r co „d is the conden- 
sation outer radius, which we take to be twice the geometrical 
mean of the deconvolved major and minor HWHM radii mea- 
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Table 4. Line-of-sight velocity dispersions and virial mass estimates for the 41 starless 1.2 mm condensations positively detected 
inN 2 H + . 



Prior to background subtraction After background subtraction 



Source 


FWHM 


Tj' 


.2mm 


cr T {p.) 


cr at 


°"AT 


M vir 




cr 




M vir 


M vir 


a- T 




a- T 






(AU X AU) 


(K) 


(M ) 


(km s - ') 


(km s _1 ) 




(M ) 




(km s~') 




(M ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


A3-MM1 


< 400 


12 


0.21 


0.21 


0.37 


1.8 


< 0.24 


< 1.2 










A-MM4 


4000 x 1400 


12 


0.29 


0.21 


0.18 


0.9 


0.60 


2.1 


0.15 


0.7 


0.52 


1.8 


A-MM5 


3700 x 2900 


12 


0.48 


0.21 


0.15 


0.7 


0.72 


1.5 


0.05 


0.2 


0.50 


1.1 


SM1N 


3000 x 1800 


20 


1.30 


0.27 


0.19 


0.7 


0.85 


0.7 


0.17 


0.6 


0.78 


0.6 


SMI 


6100 x 2100 


20 


3.20 


0.27 


0.24 


0.9 


1.57 


0.5 


0.26 


1.0 


1.68 


0.5 


A-MM6 


3200 x 2700 


20 


0.40 


0.27 


0.28 


1.1 


1.51 


3.8 










SM2 


6200 x 3400 


20 


1.30 


0.27 


0.19 


0.7 


1.67 


1.3 


0.15 


0.5 


1.43 


1.1 


A-MM8 


2900 x 2100 


20 


0.13 


0.27 


0.14 


0.5 


0.77 


5.9 


0.12 


0.4 


0.71 


5.5 


A-S 


<400 


12 


0.17 


0.21 


0.08 


0.4 


<0.07 


< 0.4 










B1-MM1 


<400 


12 


0.10 


0.21 


0.15 


0.7 


<0.09 


<0.9 










B1-MM2" 


3000 x 2100 


12 


0.17 


0.21 


0.11 


0.5 


0.46 


2.7 


0.07 


0.3 


0.40 


2.4 


B1-MM3" 


1800 x 1300 


12 


0.16 


0.21 


0.13 


0.6 


0.31 


1.9 


0.05 


0.3 


0.23 


1.5 


B1-MM4" 


4600 x 3200 


12 


0.21 


0.21 


0.13 


0.6 


0.78 


3.7 


0.06 


0.3 


0.61 


2.9 


B1B2-MM1 


2700 x 1800 


12 


0.10 


0.21 


0.15 


0.7 


0.49 


4.9 


0.07 


0.3 


0.35 


3.5 


B2-MM1 


< 400 


12 


0.14 


0.21 


0.15 


0.7 


< 0.09 


< 0.6 


0.05 


0.3 


< 0.06 


< 0.4 


B2-MM2 


4500 x 2400 


12 


0.47 


0.21 
0.21 


0.32 
0.08 


1.5 
0.4 


1.58 
0.55 


3.4 
1.2 










B2-MM3" 


< 400 


12 


0.12 


0.21 


0.22 


1.1 


< 0.12 


< 1.0 










B2-MM4 


2100 x 960 


12 


0.27 


0.21 
0.21 


0.16 
0.13 


0.8 
0.7 


0.32 
0.29 


1.2 
1.1 


0.09 


0.5 


0.25 


0.9 


B2-MM5" 


2200 x 960 


12 


0.26 


0.21 


0.20 


1.0 


0.41 


1.6 


0.12 


0.6 


0.27 


1.1 


B2-MM6 


4300 x 2700 


12 


0.78 


0.21 


0.29 


1.4 


1.47 


1.9 


0.17 


0.8 


0.82 


1.1 


B2-MM7 


<400 


12 


0.23 


0.21 


0.33 


1.6 


<0.20 


<0.9 










B2-MM8 


4000 x 4000 


12 


1.50 


0.21 


0.24 


1.2 


1.36 


0.9 


0.17 


0.8 


0.96 


0.6 


B2-MM9 


1600 x 960 


12 


0.31 


0.21 


0.27 


1.3 


0.48 


1.6 










B2-MM10 


3400 x 2200 


12 


0.60 


0.21 


0.23 


1.1 


0.90 


1.5 










B2-MM11 


< 400 


12 


0.15 


0.21 


0.43 


2.1 


< 0.31 


< 2.0 


0.26 


1.3 


< 0.15 


< 1.0 


B2-MM12 


2100 x 1300 


12 


0.39 


0.21 


0.23 


1.1 


0.52 


1.3 


0.08 


0.4 


0.28 


0.7 


B2-MM13 


<400 


12 


0.19 


0.21 


0.26 


1.3 


<0.15 


<0.8 










B2-MM14 


2100 x 1800 


12 


0.43 


0.21 


0.32 


1.6 


0.96 


2.2 










B2-MM15 


<400 


12 


0.17 


0.21 


0.14 


0.7 


<0.08 


<0.5 


0.07 


0.4 


<0.06 


<0.4 


B2-MM16 


2700 x 1300 


12 


0.35 


0.21 


0.26 


1.3 


0.70 


2.0 


0.17 


0.8 


0.45 


1.3 


B2-MM17 


<400 


12 


0.23 


0.21 


0.26 


1.2 


< 0.15 


< 0.6 


0.12 


0.6 


< 0.08 


<0.3 


C-MM2 


<400 


12 


0.12 


0.21 


0.17 


0.8 


<0.10 


<0.8 










C-MM3 


5400 x 640 


12 


0.23 


0.21 


0.12 


0.6 


0.35 


1.5 










C-MM4 


2400 x 1400 


12 


0.16 


0.21 


0.13 


0.6 


0.37 


2.3 










C-MM5 


<400 


12 


0.10 


0.21 


0.12 


0.6 


<0.08 


<0.8 


0.08 


0.4 


<0.07 


<0.7 


C-MM6 


4000 x 3700 


12 


0.33 


0.21 


0.13 


0.6 


0.76 


2.3 


0.09 


0.4 


0.65 


2.0 


C-MM7 


<400 


12 


0.13 


0.21 


0.11 


0.5 


<0.07 


<0.6 










E-MM2d 


4200 x 2700 


12 


0.63 


0.21 


0.11 


0.5 


0.61 


1.0 


0.09 


0.4 


0.58 


0.9 


E-MM4 


6900 x 5300 


12 


0.61 


0.21 


0.13 


0.6 


1.20 


2.0 


0.09 


0.4 


1.04 


1.7 


F-MM1 


4800 x 2600 


12 


0.35 


0.21 


0.19 


0.9 


0.95 


2.7 


0.13 


0.6 


0.70 


2.0 


F-MM2" 


2700 x 1600 


12 


0.17 


0.21 


0.14 


0.7 


0.43 


2.6 


0.12 


0.6 


0.39 


2.3 



(a) Condensation with two velocity components in N2H + (see TableO. Only the velocity component identified with Gaussclumps (see Table[3) 
is considered here. 

(b) For some c ondensati ons, the dust temperatures adopted here to estimate M\ lmm from the measured 1.2mm flux differ slightly from those 
assumed by MAN98. For simplici ty, we used T d = 1 2 K for all condensations except SM1N, SMI, A-MM6, SM2, A-MM8 in Oph A, which 
are likely warmer (T d = 20 K - cf. Andre et alii 19931) . The mass spectra shown in Fig.[TJand Fig.[8]were obtained using the dust temperatures 
quoted here. 



sured in the 1.2 mm dust continuum (i.e., VFWHM a x FWHMb 
- cf. col. [2] of Table S and Table 2 of iMAN98l) . For unre- 
solved condensations, the quoted M„, and a V j r values are only 
upper limits obtained under the assumption that the deconvolved 
FWHM diameter is 400 AU (i.e., ~ one fourth of the HPBW spa- 
tial resolution of the 1 .2 mm continuum observations). Note also 
that our N2H + (l-0) observations provide the line-of-sight veloc- 
ity dispersion averaged over a ~ 26" beam, while the FWHM 



angular size of the dust condensations ranges from < 15" to 
~ 30". However, since both M„, and M^mm are derived for a 
diameter ~ twice the FWHM size, our method of estimating a v ,,> 
should be reasonably accurate. Some of the secondary compo- 
nents listed in Table [2] were ignored in Table |H based on the 
Gaussclumps identifications summarized in Table [3] 

The measured N2H + (l-0) linewidths are narrow, indicative 
of small internal velocity dispersions. Prior to background sub- 
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traction, the average one-dimensional nonthermal velocity dis- 
persion, ctnt, is estimated to be ~ 0.20 ± 0.08 km s _1 among 
the 43 components listed in Table [4] More precisely, the aver- 
age one-dimensional nonthermal velocity dispersion of the con- 
densations ranges from 0.12 + 0.02 km s in Oph E, 0.13 + 
0.02 km s 1 in Oph C, 0.15 ± 0.04 km s 1 in Oph B1/B1B2, 
0.17 ± 0.04 km s 1 in Oph F, to 0.20 ± 0.08 km s 1 in Oph A, 
and 0.24 + 0.08 km s in Oph B2. For comparison, the isother- 
mal sound speed or one-dimensional thermal velocity dispersion 
is crj ~ 0.21-0.27 km s for gas temperatures of 12-20 K. 
Therefore, the condensations of Oph A, Oph Bl, Oph C, Oph E, 
and Oph F are characterized by only subsonic levels of inter- 
nal turbulence with ctnt/o't ~ 0.6-0.8 on average, while the 
condensations of Oph B2 have at most "transonic" internal tur- 
bulence with <tnt/o~t ~ 1-1 < 2 on average. This is especially 
true since the measured nonthermal component of the N 2 H + (1- 
0) linewidth may overestimate the intrinsic level of internal tur- 
bulence in the presence of infall motions and/or rotation (see 
Sect.|4~Tland|4~2l 

In contrast, the line-of-sight velocity dispersions measured 
on the scale o f the parent DCO + cores are significantly larger 
(cf. Table 2 of lLoren et al.ll 1990h . with corresponding values of 
ctnt ranging from ~ 0.17 km s _1 in Oph E, 0.2-0.3 km s _1 in 
Oph C and Oph F, to 0.35-0.5 km s" 1 in Oph A, Oph Bl, and 
Oph B2. Furthermore, the values of ctnt derived from NH3 ob- 
servations of the cores tend to b e larger, by up to a factor of 2, 
than those inferred from DCO + dLoren et alJI 1990t; A. Wootten, 
private communication). For insta nce, ctnt = 0-5 + 0.1 km s" 1 i s 
estimated in Oph Bl from NH 3 dZeng, Batrla & Wilson|[T984l) . 
compared to ctnt - 0.4 km s _1 from DCO + . Thus, the internal 
velocity dispersions derived here for the compact condensations 
are typically a factor of ~ 2 smaller than the nonthermal veloc- 
ity dispersions of their parent dense cores. We conclude that the 
level of turbulence decreases from marginally supersonic on the 
scale of the DCO + cores to transonic or marginally subsonic on 
the scale of the prestellar condensations. 

The above results are reminiscent of the transition to "co- 
herence" observed by iGoodman et alj d!998l) in isolated star- 
less cores. They suggest that, at least in nearby cluster-forming 
clouds, the initial conditions for individual protostellar collapse 
are "coherent" and largely free of turbulence. They a re also 
consis tent with scenarios, such as the "kernel" model of Mvers 
( 1998), according to which protocluster condensations form by 
dissipation of MHD tur bulence on sma ll scales within massive, 
turbulent cores (see also lNakanolll998l) . 



4. Motions within the DCO + cores of L1688 

4.1. Evidence of gravitational infall motions 

Our data s how the class i cal spectroscopic s ignature of infall mo- 
tions (cf. lEvansI 1 1999I: iMvers et ail 120001) toward at least six 
1.2 mm continuum condensations, where optically thick lines 
such as CS(2-1), CS(3-2), H 2 CO(2 12 - In), and/or HCO + (3-2) 
are double-peaked with a stronger blue peak, while low-optical- 
depth lines such as N 2 H + (101-012) and C 34 S(2-1) peak in the 
dip between the blue and red peaks of the optically thick tracers 
(cf. Fig. |5]Q. Such asymmetric line profiles skewed to the blue 



OphE— MM2d 



0phE-MM4 



0phA-SM2 



4 The optically thick spectra observed toward C-MM5 and C-MM6 
show two fainter, additional peaks at velocities ~ 2.8 and ~ 4.9 km s" 1 . 
These additional velocity components are also observed in C 18 0(2-l) 
but are not seen in N9H + (l-0). They likely arise from the ambient cloud 
rather than from C-MM5 et C-MM6 themselves. 
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Fig. 5. Spectroscopic signatures of infall motions observed to- 
ward 6 Ophiuchus starless condensations: (a) OphE-MM2d, (b) 
OphE-MM4, (c) OphA-SM2, (d) OphB2-MM16, (e) OphC- 
MM5, and (f) OphC-MM6. On each panel, the vertical dot- 
ted line marks the source systemic velocity as derived from a 
Gaussian HFS fit to the observed N2H + (l-0) spectrum. 



in optically thick tracers are produced when a gradient in excita- 
tion temperatu re toward source cent er is combined with inward 
motions (e.g. lLeung & BrowrJI 19771) . Here, among a total of 25 
starless condensations observed in at least one optically thick in- 
fall tracer and one optically thin transition, blue infall profiles 
are clearly observed in 6 condensations (A-SM2, B2-MM16, C- 
MM5, C-MM6, E-MM2d, and E-MM4, see Fig. g), tentatively 
observed in 10 other condensations (A-MM4, A-MM5, A-SM1, 
A-MM8, C-N, B2-MM12, B2-MM13, B2-MM14, B2-MM15, 
B2-MM17), and not seen in the remaining 9 sources (A-MM6, 
A-SM1N, F-MM1, F-MM2, B1-MM4, B2-MM2, B2-MM6, B2- 
MM8, B2-MM10). 

In the case of A-SM2 and C-MM6, the absorption dip ob- 
served in the optically thick line profiles is slightly redshifted 
with respect to the LSR velocity traced by the optically thin 
lines, by ~ 0.1 and ~ 0.1-0.2 km s _1 , respectively. This sug- 
gests that the outer gas layers responsible for the absorption dip 
in these two condensations are characterized by significant in- 
ward velocities ~ 0.1 km s _1 (see, e.g., discussion in Sect. 3.3 of 
Bell oche et al.ll 2002). By contrast, the absorption dip is not red- 
shifted (nor blueshifted) toward E-MM2d, E-MM4, B2-MM16, 
and C-MM5, suggesting that only comparatively deeper layers 
undergo infall motions in the latter sources. 
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Fig. 6. Centroid velocity maps of Oph B2 (a), Oph Bl (b), Oph A (c), Oph F (d), Oph E (e), and Oph C-S (f) as derived from 
HFS fits with one or two velocity components to the N2H + (l-0) spectra with a signal-to-noise ratio larger than ~ 5 (filled squares 
of varying sizes and colors). The (0,0) positions are given in the caption of Fig. [2] The linear size of the squares increases, and their 
color changes from blue to red, as V/ sr increases. The underlying contours represent the same N2H + (l-0) integrated intensity maps 
as shown in Fig.|2^-f. The direction of the fitted velocity gradients in Oph F, Oph C-S, and Oph E are shown by arrows in panels 
(d), (e), and (f) (see Table[5]). In Oph F, two velocity components were fitted separately and the arrows show the two corresponding 
velocity gradients associated with components Oph Fl at v ~ 4.2 km s -1 and Oph F2 at v ~ 4.6 km s _1 , respectively. Crosses 
mark the 1 .2 mm continuum positions of starless condensations (MAN98), while stars mark the positions of protostars (VLA 1623, 
YLW 15/IRS 43, CRBR 85, and WL 19). 



4.2. Large-scale velocity structure 

In order to investigate the presence of systematic velocity gra- 
dients, such as rotational gradients, within the DCO + cores of 
L1688, we plot maps of the N2H + (l-0) centroid velocity across 
Oph A, Oph Bl, Oph B2, Oph C-S, Oph E, and Oph F in Fig.gh- 
f. These centroid velocity maps were derived from the origi- 
nal N2H + (l-0) data cubes before background subtraction. No 
clear, large-scale velocity gradient can be seen in the velocity 
maps of Oph A, Oph Bl, and Oph B2, and the underlying kine- 
matic structure cannot be described as a simple superposition 
of a small number of velocity gradients. In Oph A, which re- 
sembles a curved filament in the integrated intensity map, the 
two ends of the filament are characterized by lower systemic 
velocities (V/ s , ~ 3.2 and ~ 3.5 km s _1 ) than the central re- 
gion (V/ sr ~ 3.7 km s _1 ), dominated by the objects A-SM1 and 
VLA 1623. In Oph B2, the two condensations B2-MM10 and 
B2-MM15 stand out with significantly higher systemic veloci- 
ties (V/ sr ~ 4.4 km s~') than the bulk of the core (Vi sr ~ 3.7-4.1 
km s~'). In Oph Bl, two velocity components are detected on 
the lines of sight to B1-MM2, B 1 -MM3, and B 1 -MM4. Based on 
the results of our Gaussclumps analysis (see Tabled, we believe 
that only the ~ 4.0 km s component, connected to Oph B2, is 
physically associated with the compact condensations B1-MM2, 
B1-MM3, andBl-MM4. 



Well-defined N2H + (l-0) velocity gradients were identified 
toward Oph C-S, Oph E, and Oph F, and their amplitudes es- 
timate d using the least-squares fitting method of iGoodman et alj 
dl993l) . In particular, Oph C-S is characterized by a system- 
atic projected velocity gradient of mean value ||VV|| ~ 4.4 
km s -1 pc , with velocities increasing along the direction of 

position angle P.A. 77° in the plane of the sky (see Table [5] 

and Fig.|6f). 

In Oph E, our N2H + (l-0) data suggest the presence of a ve- 
locity gradient ||VV|| ~ 3.0 km s" 1 pc -1 at P.A. 101° (cf. 

Fig-B)- The signal-to-noise ratio in our N2H + (l-0) map is how- 
ever insufficient to constrain the nature of the velocity field. 
Fortunately, we can also use our higher signal-to-noise maps 
taken in H I3 CO + (1-0), DCO + (2-l), and DCO + (3-2) to further 
constrain the velocity structure. For Oph E as a whole, the three 
tracers indicate a mean velocity gradient 5.2 + 0.3 km s _1 pc~' 
at P.A. = -108° ± 4°, in fairly good agreement with the gradient 
suggested by the N 2 H + (l-0) data (cf. Fig.|6t). 

In Oph F, our N2H + (l-0) data indicate the presence of two 
distinct velocity components, overlapping near the position of 
F-MM2 in the plane of the sky (see Fig. |6}i). Oph F may thus 
actually consist of two independent dense cores, Oph Fl and 
Oph F2, with differing line-of-sight velocities (~ 4.2 and ~ 4.6 
km s~'). Figure[6]i shows the direction of the mean velocity gra- 
dient derived in each of these two cores (see Table [5] for derived 
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Size" 
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II v v || 


P.A. d 








(" x ") 


(km s" 1 ) 


(km s'pc~^) 


o 


(km s" 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Oph C-S 


N 2 H + (l-0) 


120 x 70 


14 


3.72 


4.4 


-77 


0.042 


OphFl 


N 2 H + (l-0) 


150x50 


16 


4.17 


3.8 


17 


0.026 


OphF2 


N 2 H + (l-0) 


50 x 80 


5 


4.60 


6.3 


-119 


0.028 


OphE 


H 1J CO+(1-0) 


180 x 80 


19 


4.47 


5.0 


-106 


0.052 




DCO + (2-l) 


170 x 70 


43 


4.53 


5.1 


-104 


0.051 




DCO + (3-2) 


160 x 60 


68 


4.48 


5.6 


-113 


0.044 




Average 








5.2(3) 


-108(4) 





Notes: the numbers in parentheses indicate the uncertainty in units of the last digit. 

(a) Approximate diameter of the region over which a velocity gradient was fitted 

(b) Number of fully independent points used in the gradient fitting. 

(c) The centroid velocity maps were fitted with the function Vo + WAX, with AX the position vector measured from the (0,0) position. 

(d) Position angle of the direction of the velocity gradient W. 

(e) rms residual of the fit. 



parameters). The fact that the velocity gradient found in Oph Fl 
(~ 3.8 km s pc~' at P.A. ~ 17°) is nearly opposite to the gra- 
dient in Oph F2 (~ 6.3 km s" 1 pc" 1 at P.A. ~ -1 19°)0 supports 
the view that Oph F is made up of two independent cores, partly 
overlapping in projection, and not a single core with a mean ve- 
locit y gradient ~ 1.8 km s -1 pc~' at P.A. ~ -50° as suggested 
bv lLoren etaT] d!990l) based on lower-resolution DCO + (2-l) ob- 
servations. 

For L1688 considered as a whole (cf. Fig. [7J, a global ve- 
locity gradient of ~ 1.1 km s _I pc~' is seen from North- West 
(Oph A) to South-East (Oph F) (see also Loren 1989 in 13 CO 
and Loren et al. 1990 in DCO + ). The direction of this large- 
scale gradient is P.A. ~ 120° (cf. note a of Table [6j, as deter- 
mined by a least-squares fit to the observed distribution of the 
line-of-sight systemic velocities of the condensations detected 
in N2H + (l-0). The observed velocity pattern is however clearly 
more complex than that exp ected from simple large-scale rota- 
tion (see also lLoren|fl989l) . For instance, the direction of the 
mean velocity gradient changes to P.A. ~ 180° if the condensa- 
tions of Oph Bl and Oph B2 are ignored. 

4.3. Relative motions of the protocluster condensations 
within the DCO + cores 

Our estimates of the centroid line-of-sight velocities toward the 
prestellar condensations of L1688 (cf. col. 7 of Table[2j provide 
interesting constraints on the relative motions between conden- 
sations. Figure [7] shows the distribution of systemic velocities 
as derived from our N2H + (l-0) HFS fits, overlaid o n the lowest 
contours of the 1 .2 mm continuum map of lMAN98l 

Both the velocity differences between neighboring conden- 
sations and the overall velocity dispersion of the condensations 
within the L1688 protocluster are small. Table [6] gives the one- 
dimensional velocity dispersion, <Tw, derived from our N2H + (1- 
0) observations for various samples of condensations (cf. col. [6] 
and col. [7]). This velocity dispersion was estimated in two 
ways, first around the mean LSR velocity of the condensations 
(estimate given in col. [6]), and second around the mean centroid 
N2H + (101-012) velocity measured in the maps of Fig. [6] (esti- 
mate given in col. [7]). For the ensemble of 41 compact prestel- 

5 The direction of the velocity gradient is however more uncertain in 
Oph F2 since it was derived using only 5 independent points. 



lar condensations and 3 protostars detected in N2FL, we estimate 
a global velocity dispersion cr 1D ~ 0.36 km s _1 about the mean 
systemic velocity <V/ Jf >~ 3.95 km s _1 . Assuming isotropic ran- 
dom motions, this corresponds to a three-dimensional velocity 
dispersion cr 3z > = V3 ctw ~ 0.62 km s 1 (cf. col. [8] and 
[9] of Table |6). If the velocity distribution is Maxwellian, the 
mean condensation speed relative to the center of mass of the 
system is V meml = V8/tt cr ID ~ 0.57 km s _1 and the mean 
relative speed between condensations is V re i - V2 V mea „ - 
<T\ D ~ 0.81 km s Adopting a diameter D ~ 1.1 pc 
re.g. lWilking&Ladall983h for the L1688 cloud, such a velocity 
dispersion implies a typical crossing time, f cr0JJ = Djcr^ ~ 1.8 
xlO 6 yr (cf. col. [10] and col. [11] of Tabled), for the condensa- 
tions within the L1688 protocluster. After subtracting the global, 
systematic velocity gradient of ~ 1.1 km s _I pc~' seen across 
L1688 (cf. Sect. |4.2| above), the resulting ID velocity dispersion 
of the condensations is even lower, <x IO ~ 0.25 km s _1 , suggest- 
ing that the crossing time associated with purely random conden- 
sation motions is larger, t cross ~ 2.5 x 10 6 yr. The crossing times 
estimated separately within the individual DCO + cores with a 
statistically significant number of condensations (i.e., Oph A, 
B2, C) are only a factor of ~ 2 shorter, t cross ~ 0.6-1.0 xlO 6 yr. 

5. Discussion 

5.1. Dynamical state and fate of the L1688 condensations 

The starless condensations identified by MAN98 in L1688 are 
highly centrally concentrated and feature large density contrasts 
over the local background medium, strongly suggesting they are 
self-gravitating. More specifically, the estimated mean densities 
of the condensations exceed the mean densities of the parent 
DCO + cores by a typical factor ~ 5-20, while the mean col- 
umn densities of the condensations exceed the background core 
column densities by a factor i 2. For comparison, a critical 
self-gravitating Bonnor-Ebert isothe rmal spheroid has a mean 
density contrast pBE/Pext ~ 2.4 (e.g. lLombardi & Bertinll200Th 
and a mean column density contrast £BE/2 ex t ~ over tne 
external medium. 

The line observations reported in this paper provide more 
direct evidence that most of the L1688 condensations are grav- 
itationally bound. Indeed, the narrow linewidths measured in 
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Fig. 7. Line-of-sight systemic velocities of the 41 preste llar condensations L1688 detected in N2H + (l-0), overlaid on the lowest 
contours of the 1.2 mm continuum mosaic of MAN98. These velocities correspond to those listed in col. [7] of Table [2] Each 
condensation is represented by a filled circle whose size increases with Vlsr. The color coding varies from black to yellow with 
increasing Doppler shift. The systemic velocities of the protostars VLA 1623, YLW 15/IRS 43 and CRBR 85 are also shown. A few 
condensations have two markers, reflecting the presence of two velocity components in their spectra. 



N2H + (l-0) (see Sect. 13.2b imply virial masses which generally 
agree wi thin a factor of ~ 2 with the mass estimates derived by 
MAN98 from the 1 .2 mm dust continuum. This can be consid- 
ered a good agreement since both M V! > and M\2mm are them- 
selves uncertain by a factor of ~ 2. On the theoretical side, self- 
gravitating condensations are expected to have virial mass ra- 
tios a V i r = M V i r jM\2mm ~ 2, while obj ects in gravitational viria l 
equilibrium should have a n > ~ 1 (e.g. lBertoldi & McKedl 1992k 
Here, among the 43 velocity components measured toward the 
41 compact condensations positively detected in N2H + (l-0), 15 
have an estimated virial mass ratio a V! > < 1 and another 23 have 
a V i r i 2.5. Altogether, 38 components have ar V! > i 2.5 and only 5 
have a V i r i 3. Among the latter 5 components, only 2 are associ- 
ated with condensations more massive than M\2mm ~ 0.35 M , 
of which one is associated with a condensation (B2-MM2) hav- 
ing another velocity component with a Ylr < 2. Therefore, 37 of 
the 41 compact 1.2 mm condensations detected in N2H + (l-0) 
have at least one velocity component with a v / r t 2.5. 

We conclude that a large majority (~ 75-90%) of the L1688 
condensations detected in N2H + (l-0), including ~ 90% of those 
more massive than ~ 0.35 M , are likely gravitationally bouncQ. 

6 Note that the surface pressure term of the virial theorem, which we 
neglect here, tends to reduce the mass required for virial equilibrium, 
i.e., to increase the value of a vir expected in virial equilibrium - see 
lBertoldi&McKeeHl992h . 



Coupled with the detection of signatures of infall motions toward 
some of them (see Sect. 14. IK this support s the noti on that ~ 75% 
of the starless condensations identified by |MAN98l at 1 .2 mm are 
prestellar in nature a nd on the verge of form ing protostars (cf. 
lAndre et al.1 2000 and IWard-Thompson et all 2007 for a defini- 
tion of the prestellar stage). A fraction (~ 30%) of the conden- 
sations with Mi 2mm ~ 0.35 M Q may be only marginally bound. 
The status of the condensations less massive than ~ 0. 1 M is 
less clear since these are weaker and often undetected in N2FL. 
Some of them may possibly correspond to unb ound transient 
objec ts generated by supersonic turbulence (cf. iKlessen et al.l 
2005). Interesti ngly, if a significant fraction of the condensations 
identified by |M AN98I below ~ 0.35 M Q are indeed unbound and 
not truly prestellar in nature, then the apparent excess of star- 
less condensations at the low-mass end of the mass spectrum 
compared to the IMF of field systems (cf. Fig. 1) may find a nat- 
ural explanation in the context of a one-to-one mapping between 
prestellar condensations and stellar systems. 

5.2. Lifetime of the L1688 condensations 

The lifetime of the L1688 prestellar condensations is rather un- 
certain but several lines of reasoning suggest a typical value of 
~ 10 5 yr. Table [7] provides estimates of some relevant evolu- 
tionary timescales. First, based on the measured mean densi- 
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Table 6. Velocity dispersion of the L1688 protocluster condensations. 
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Sample" 


n b 


D c 


<vf> 

Isr 




T 


<d,o 


< 


3D,c 


D/o-3 D 


D/0-3D.C 






(pc) 


(km s _1 ) 


(kms -1 ) 


(kms^ 1 ) 


(kms" 1 ) 


(kms 4 ) 


(km s" 1 ) 


(10 s yr) 


(10 6 yr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Oph A 


9 


0.28 


3.44 


3.48 


0.19(5) 


0.20(5) 


0.33(8) 


0.34(9) 


0.8 


0.8 


OphBl/B2 
Oph C,E,F 


24 
10 


0.33 
0.44 


4.05 
4.09 


3.96 
4.05 


0.20(3) 
0.39(9) 


0.22(3) 
0.39(9) 


0.35(5) 
0.67(16) 


0.38(6) 
0.68(16) 


1.0 
0.7 


0.9 

0.6 


L1688 


47 


1.10 


3.95 


3.90 


0.36(4) 


0.36(4) 


0.62(7) 


0.63(7) 


1.8 


1.7 


L1688 -VV 


47 


1.10 






0.25(3) 




0.43(4) 




2.6 





Notes: the numbers in parentheses indicate the uncertainty in units of the last digit. 

(a) The first three samples contain the velocity components of Table|4]only. The last two samples include the protostars in addition. For the 
"L1688 -VV sample, the velocity dispersion is computed after removing the large scale velocity gradient measured with the method described 
in note cof Tablegjusing C-N as the (0,0) position: V Q = 3.87 km s" 1 , ||VV|| = 1.1 km s^pc" 1 , P.A. = 117 °. 

(b) Number of velocity components used for the calculations. 

(c) Diameter of the region containing each sample. 

(d) Mean LSR velocity of the components, computed with the velocities listed in col. [5] of Table|2] 

(e) Mean centroid velocity computed on the N2H + ( 101-012) spectra shown in Fig. [6] 

(/) Standard deviation of the distribution of component LSR velocities around <V/ Jr > given in col. [4]. The error bar was estimated as 

o"id/ V2(n - 1), assuming that the sample is drawn from a larger population whose velocity distribution follows Gaussian statistics, 
(g) Same as in col. [6] but computed around <V cen f> given in col. [5]. 

(/j) 3D velocity dispersion calculated from cr m assuming isotropic motions. The error bar was scaled from that estimated for <x 1D . 
(0 Same as in col. [8] but computed around <V ce „,> given in col. [5]. 



ties «H2 = pll*tn H ~ 10 s - 10 7 cm 3 , the free-fall dynamical 
timescales of the condensations, tff = (3n/32Gp) 1 ^ 2 , range from 
t ff ~ 10 4 yr to tff ~ 10 5 yr. While some of the condensations 
show evidence of infall motions, most of them are unlikely to be 
in free-fall collapse. We may conservatively estimate their evo- 
lutionary timescale to be ~ 3 trt, which is appropriate for mag- 
netica lly supercritical dense cores (e.g. ICiolek & Mouschoviasl 
1 1994b and typical of is olated prestellar c ores detected in the sub- 
millimeter continuum dKirk et al.l l2005). This leads to timescale 
values ranging from ~ 2 x 10 4 yr to ~ 3.5 x 10 5 yr (cf. col. [6] of 
Tabled- Second, assuming that the L1688 condensations evolve 
into stellar systems at a constant rate, we may derive a rough 
statistical estimate of their lifetime (col. [7] of Tabled by com- 
paring the number of observed prestellar condensations (col. [3] 
of Table 13 to the number of pre-main sequence (PMS) systems 
found in the sa me region (col. [21 of T able Q. The ISOCAM 
mid-IR survey of Bontemps et al. (2001) revealed a total of 109 
Class II PMS objects in L1688 with ages ~ 0. 4 - 1 Myr . If all 
of the 45 starless condensations identified by IMAN98I above 
0.1 M Q are truly prestellar in nature, this points to a prestel- 
lar lifetime ~ 2 - 4 x 10 5 yr. A third timescale estimate may 
be obtained by dividing the typical condensation outer radius 
r C on d ~ 4000 AU by a typical infall speed V- m f ~ 0.1 - 0.3 km/s 
(cf. lBelloche. A ndre & Motte 2001). This approach yields a con- 
densation lifetime ~ 0.6-2 xlO 5 yr. We conclude that the 1.2mm 
continuum condensations of L1688 are likely characterized by a 
range of lifetimes between ~ 2 x 10 4 yr and ~ 5 x 10 5 yr (see 

col. [6] and col. [7] of Table |7). 

As pointed out by IClark etail d2007l) . if the lifetime of the 
condensations depends on their mass, then the observed mass 
spectrum is not necessarily representative o f the intrinsic con - 
densation mass distribution (CMD) (see also Elmegreen 2000). 
This is due to the fact that an observer is more likely to detect 
long-lived condensations than short-lived condensations. Here, 
however, the mean densities of the L1688 condensations are es- 
sentially uncorrected with their masses, so that there is no sys- 
tematic dependence of the dynamical timescale on the mass. To 
quantify the importance of the potential timescale bias, we plot, 
in Fig. [8j a weighted version of the central Ophiuchus CMD 




Mass, M (M ) 



Fig. 8. Weighted cumulati ve mass s pectrum of the 57 starless 
condensations identified bv lMAN98l (histo gram with error bars), 
compared to the same mass distributions as shown in Fig. Q] as 
well as the Salpeter power-law IMF (solid line). Here, each con- 
densation was assigned a weight inversely proportional to its es- 
timated free-fall dynamical timescale (see text). The flattening 
apparent below ~ 0.4 M in the unweighted CMD (Fig. [TJ is 
not seen in the weighted CMD shown here, which is essentially 
consistent with a single, Salpeter-like power law. 



in which each condensation was assigned a weight equal to 

~ifs~ = <% 2 > ( mstea d °f 1 as use d f° r Fig. [U, where < tff > is 
the average free-fall time of the condensations. Such a weight- 
ing should allow us to recover the intrinsic shape of the CMD 
assuming that the lifetime of each condensation is proportional 
to its free-fall time. As can be seen in Fig. [8] this weighting does 
not change the high-mass end of the CMD and only affects the 
low-mass end. We conclude that the steep, Salpeter-like slope of 
the CMD at the high-mass end is robust, but that the presence of 
a break at ~ 0.4 M is less robust. 
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Table 7. Evolutionary timescales for various systems of protocluster condensations. 
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5.3. Likelihood of interactions between condensations 

Using the velocity dispersions derived for the condensations in 
Sect. 14.31 and the lifetime estimates of Sect. 15.21 above, we are 
now in a position to assess whether dynamical interactions be- 
tween condensations are a frequent or rare phenomenon in the 
L1688 protocluster. 

Assuming an isotropic, Maxwellian velocity distribution, the 
time, t co u, required for a condensation to interact or collide with 
another condensation can be estimated from the following for- 
mula for the collision rate l/t co u : 

I /ten = 4 V^r x n cond cr 1D r 2 cond x (1 + 0) , (1) 

where n com i is the number density of condensations, r con d is 
the condensation outer radius (which we take to be twice the 
measured HWHM radius at 1.2mm - cf. Sect. 13.21 ). nr 2 , is 
the collisional cross section of the condensation, and 1 + = 
1 + GM con dl{o\ D r C ond) is the gravitational focusing factor (e.g. 
Binn ev & Tremairielll987l) . Further assuming that the conden- 
sations are homogeneously distributed in a spherical system of 
radius R, so that n con d - N CO nd/(jnR 3 ), the ratio of t co n to t cross 
takes on the simple form: 

t co n 1 fjt R 2 1 

= - A /-X — X . (2) 

t cross 2 V 3 Ncond r- cond 1 + © 

Using the above equation for the entire L1688 cluster- 
forming clump (R ~ 0.55 pc, N CO nd = 57, cr 1D = 0.36 km s _1 ) 
and adopting typical condensation properties (r con d ~ 2500 AU 
and Mcond ~ 0.4 M ), we find 8 /^ 88 ~ 9 and thus t£f 8 ~ 
16 X 10 6 yr. The latter is two orders of magnitude larger than the 
estimated lifetime of the condensations (Sect. 15721 and about one 
ord er of magnitude longer than the age of the L1688 PMS clus- 
ter dBontemps et al.ll2001l) . For simplicity, our analysis ignores 
the dynamical impact of the ambient protocluster environment. 
In such a dense environment, gas drag is likely to have an effect 
on the motions of individual condensations, which will make the 
associated velocity field depart from a true Maxwellian distribu- 
tion. A detailed assessment of the effects of gas drag is beyond 
the scope of this paper, but qualitatively at least we expect the 
environment to make the interaction timescales even longer than 
the above estimates. We conclude that the interaction process is 
much too slow on the scale of the entire L1688 system to play 
any significant role in the evolution of the protocluster conden- 
sations. 

Dynamical interactions are however more likely to occur on 
smaller scales than the whole L1688 cl uster-forming clump. In 
their mid-IR census of the PMS cluster, Bontemps et Id] d2001l) 
identified three main sub-clusters associated with the DCO + 
cores Oph A, Oph B1/B2, and Oph C/E/F, respectively. The fact 
that sub-clustering is observed in the spatial distribution of PMS 



objects further supports our conclusion that significant interac- 
tions between individual objects have not yet occurred on the 
scale of the entire protocluster. Using Eq. (2) for the conden- 
sations associated with each of the three sub-cluster systems 
(R ~ 0.15 pc, N coni i = 10 — 20), we find shorter interaction 
timescales than for L1688 as a whole, f™* ~ 6 - 22 x 10 5 yr 
(see col. [9] of Table |7). The derived collision timescales nev- 
ertheless remain significantly longer than the typical condensa- 
tion lifetime ~ 10 5 yr, indicating that dynamical interactions 
between condensations cannot be a dominant process even on 
the scale of the sub-clusters and DCO + cloud cores. The colli- 
sion timescales are even longer than the combined lifetime of 
the prestellar and protostellar (Class 0/Cl ass I) phases, which we 
estimate to be at most ~ 5 x 10 5 yr (e.g. lGreene et al.ll 1 9941) . In 
general, therefore, the prestellar condensations do not have time 
to orbit through their parent DCO + core and collide with one 
another before evolving into Class II PMS objects. We note that 
tcoll is only slightly larger than 5 x 10 s yr in the case of the Oph A 
cluster-forming core, suggesting that a few interactions may oc- 
cur during the whole evolution of this sub-cluster. Interestingly, 
based on a recent detailed study of the N GC 2264-C cluster- 
forming clump in the Mon OBI complex, Peretto et al. (2006) 
found direct evidence of a strong dynamical interaction in the 
dense inner region of that protocluster. However, the interaction 
observed at the center of NGC 2264-C is purely gravitational in 
origin and results from large-scale, co herent collapse moti ons as 
opposed to random turbulent motions ( Peretto et al.l l2007). 

Our present results in L1688 seem inconsistent with models 
which resort to strong dynamical interac tions to build up a mass 
spectr um comparable to the IMF (e.g. iPrice & Podsiadlowskil 
1995 - see also Sect.l574lbelow). 

5.4. Comparison with the competitive accretion picture 

In a scenario of clustered st ar formation propose d by Bonnell, 
Bate and collaborators (e.g. iBonnell et al.ll200T7allbl) . the IMF is 
primarily determined by competitive accretion and dynamical 
interactions/ejections during the protostellar phase, correspond- 
ing observationally to Class and Class I objects. In this pic- 
ture, turbulence generates density fluctuations within molecu- 
lar clouds, some of which are gravitationally unstable and col- 
lapse to interacting protostars or protostellar seeds. The initial 
envelope/core masses of these protostars are essentially uncor- 
rected with the corresponding fin al stellar masses, especially a t 
the high-mass end of the IMF (cf. lBonnell. Vine. & Batd l2004). 
The protostars acquire most of their mass by moving around 
in the gravitational potential well of the parent cluster-forming 
clump and accreting background gas that initially did not be- 
long to the corresponding protostellar envelopes/cores. This pro- 
cess of competitive accretion of background gas is highly non- 
uniform and depends primarily on the initial stellar position /?* 
within the protocluster. The few protostars initially located near 
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Table 8. Rates of competitive mass accretion for a range of objects in the L1688 protocluster. 
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3500 
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3000 


6 x 10~ 7 


Class in inner Oph A core 


4 x 10 5 


0.3 


2500 
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Class from multiple system 
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0.2 


2000 


3.5 x 10~ 6 


inside collapsing condensation 










Condensation in inner Oph A core 


4x 10 5 


0.3 


2500 


3 x 10~ 6 



the center of the cluster potential accrete rapidly from the start 
and become massive stars, while protostars in the low-density 
outer regi ons accrete much m ore slowly and become low-mass 
stars (see iBonnell et al.ll2001bl for a quantitative toy model con- 
sistent with the observed IMF). 

In the context of this scenario, it is difficult to explain why 
the Ophiuchus prestellar condensations have a mass distribution 
resembling the IMF. Our present results on the kinematics of the 
protocluster condensations allow us to further discuss the possi- 
ble relevance of the competitive accretion picture in the case of 
the L1688 embedded cluster. The rate at which a protostar ac- 
cretes mass competitively as it travels through th e background 
protoc luster gas is M acc * npi, ack v re i R 2 acc (cf. IBonnell et alj 
2001a), which gives: 



Macc ~ 1.2x10 M yr" 



-back 
'H2 



10 5 . 



Vrel 



0.3km/s 



/ Race \ 2 

\3000AU/ ' 



(3) 



where pback is the background gas density, v re i is the relative 
velocity between the protos tar and the lo c al am bient gas, and 
Race is the accretion radius. IBonnell et all d200 lab have shown 
that a good analytic approximation for R acc is provided by the 
smaller of the Bondi-Hoyle radius, Rbh = 2GM+/ (v 2 rd + cf), 

(M \l/3 
XfM R*, where 

c s is the gas sound speed and M enc (/?*) is the mass enclosed 
within the protocluster at the protostar's position The Bondi- 
Hoyle radius, Rbh, is the radius where the gravitational po- 
tential due to the protostar exceeds the kinetic energy of the 
gas. Here, we estimate the relative gas-protostar velocity to be 
the derived mean condensation speed V mea „ = cr\ D ~ 

0.3 - 0.4 km s _1 < 2c s , assuming a Maxwellian velocity dis- 
tribution. Therefore, the typical value of the Bondi-Hoyle radius 
is: R BH ~ 9900 AU (M+/0.5 M ) (v re ;/0.3 km/s)" 2 . 
The tidal radius, R t idah expresses the fact that the tidal forces 
exerted by the gravitational potential of the ambient protocluster 
limit the zone of influence of a given protostar. In gas-dominated 
protoclusters such as L1688, one usually has Rtidal ~ Rbh and 
thus R acc 

Rbh ~ Rtidal and thus R acc * R BH (IBonnell et al.ll2001al) 

To estimate the typical value of R t uiai in the central 
Ophiuchus case, we assume that the overall gas distribution in 
the centrally-condensed L1688 protocluster and individual sub- 
clusters (e.g. Oph A) follows a p oc r~ 2 density profile on av- 
erage. Thi s assump tion is consistent with the density gradient 
derived by MAN98 for the outer parts of the DCO + cores based 
on the circularly-averaged profiles observed at 1.2mm. Under 
this assumption, M mc has a very simple expression, M enc (R+) = 

M c i us X(RJRclus), and thus R tidal * 0.5 (^zf 3 «J£ R 2 J\ where 
M c ius and R c [ us are the total (gas + stars) mass and outer radius of 
the protocluster, respectively. For the entire L1688 protocluster, 



we have M dus ~ 650 M Q , R c i us ~ 0.55 pc, and thus ~ 

3500 AU (M*/0.5M G ) 1/3 (fl*/0.3pc) 2/3 . For the Oph A sub- 
cluster, we adopt M c i us ~ 30 M , R c i us ~ 0.14 pc, and thus 
Kdai ~ 3000 AU (M*/0.5M o ) 1/3 (.R*/0.1pc) 2/3 . As expected, 
we see that R tida i < R BH and thus R acc as R lida i ~ 3000 AU. 



Based on these simple estimates, we find that the typical 
tidal-lobe radius of the L1688 condensations/protostars is com- 
parable to the actual radiu s of the c ondensations and protostellar 
envelopes as measured by MAN98 at 1 .2 mm. This suggests that 
the overall gravitational potential of the L1688 protocluster does 
play an important role in limiting the condensation/envelope 
masses, in qualitative agreement with the competitive accretio n 
picture (cf. IBonnell eTall 12001 at IBonnell. Vine. & Batdl2004l) . 
Note that the small value found here for the mean condensa- 
tion velocity relative to the local gas (V mea „ ~ 0.3 km s ) is 
also consistent with the predictions of the competitive accretion 
model during the gas- dominated phase of protocluster evolution 
(IBonnell et al.ll2001al) . However, when we estimate the mass ac- 
cretion rate resulting from competitive accretion in various typi- 
cal situations, we find relatively low values (see Table [3), which 
are smaller than the infall rate expected from the gravitational 



collapse of individual conde nsations, i .e., M,„/ ~ 1 — 10 c s /G 
(see IShul 1977 and iFoster & Chevalier! 19 93 for model predic- 
tions, and iBelloche. Hennebelle & Andrei 2006 for an observed 
example). For a Class protostar such as VLA 1623 embedded 
in the inner part of the Oph A DCO + core, we estimate the back- 
ground gas density to be relatively high, njg* 
and the tidal radius to be at most /? v/ r A , 1623 

tidal 

3 x 10~ 6 



'M yr 



2500 AU, giving 
For comparison, the infall rate due to 



gravitational collapse is significantly hi gher at this early stage 
M inf > 10 c 3 /G > 10~ 5 - 10" 4 M yr" 1 feontemps et all 119961: 
iJavawardhana et al.ll200lt lAndre. Motte & Bellochell200ll) . For 
a more evolved Class I protostar such as YLW 15/IRS 43 em- 
bedded in the Oph F DCO + core, the background gas den- 



sity is lower, n 
at most R 1R f^ 

tidal 



back 
H2 



5 x 10 cm , and the tidal radius is 



„ , 3000 AU, giving M acc ~ 6 x lO- 7 M yr- 
r V ! rS D ely ' dUSterS ** comparison, the gravitational infall rate is M inf ~ c 3 /G 



2 x 10 6 M yr 1 at this late protostellar stage (e.g. Adam set al.l 
119871: iBontemps etailll996l) . For both Class and Class I ob- 
jects we thus find M acc < M{ n f/3, implying that local gravita- 
tional collapse domin ates over competitive accretion (see also 
iRrumholz et alj|2005l) . We conclude that competitive accretion 
at the protostellar stage is unlikely to be the main mechanism re- 
sponsible for determining the final masses of stellar systems in 
the central Ophiuchus star-forming cloud. If each prestellar con- 
densation frag ments into a small-N system during protostellar 
collapse (e.g. iGoodwin et alj|2007l) . then competitive accretion 
may possibly play a more important role in defining the final 
masses of the individual components. 
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We also believe that competitive, Bondi-like accretion is 
more likely to operate at the prestellar stage and may pos- 
sibly govern the growth of starless condensat i ons within a 
cluster-forming clo ud (cf . iBonnell et al.l 1200 lbt iMversi l2000t 
Basu & Jones 2004). The growth rate of a typical ~ 0.5 M con- 
densation embedded in the inner part of a DCO + core such as 
Oph A should be similar to the competitive accretion rate esti- 
mated above for VLA 1623, i.e., M acc ~ 3 x 10~ 6 M yr _1 . This 
is sufficient to approximately double the mass of the condensa- 
tion in ~ 2 x 10 5 yr, a timescale comparable to the condensation 
lifetime. However, once fast, nearly free-fall collapse sets in, we 
expect the growth rate to be quickly overwhelmed by the infall 
rate, so that a collapsing prestellar condensation should not have 
time to grow significantly in mass before evolving into a PMS 
system. 

Another feature of the dynamics of a gas-dominated proto- 
cluster is that the entire system is expected to undergo global 
collapse/contraction, resulting in a centrally-condensed over- 
all structure much like a self-grav itating isothermal sphere (cf. 
lAdamsll2000l IBonnell et all 12001 bl) . In this view, both the gas 
accretion and the protocluste r evolution occur on the global dy- 
namical timescale (see also Klessen & Burkert 2000l). 
Observationally, the velocity dispersion, <x 1D , estimated in 
Sect. 14.31 above for the L1688 condensations implies a bind- 
ing virial mass M,,,> « 3 x Rcr 2 w IG ~ 50 M Q , which is 
much less than the total gas mass ~ 550 M H of the associ- 
ated C ls O cloud (cf. lWilking&Ladalll983h and even less than 
the total stellar mass of the present infrared embedded clus- 
ter (M stars ~ 100 Mm. iBontemps et ail 1200 lb . In other words, 
the observed condensation-to-condensation velocity dispersion 
is a factor of ~ 3 smaller than that expected in virial equilib- 
rium (see also Table 6 of iPeretto et al.l 120061) . Since the mag- 
netic field does not seem to b e strong enough to support the 
cloud (e.g. iTroland et al.lll996l) . this comparison suggests that 
the L1688 system is indeed gravitationally unstable and possi- 
bly in an early state of large-scale, magnetically supercritical 
contraction. Interestingly, the CO and 13 CO lines observed to- 
ward the L1688 cloud exhibit the classical spectroscopic signa- 
ture of contraction motions (cf. Sect. |4.1| above) over most of the 
protocluster extent (see lEncrenaz et alll 19751) . We also note that 
L1688 is not the only protocluster for which both subvirial rel- 
ative speeds and evidence of large-scale co ntraction have been 
found . Two other exa mples are NGC 2264 (IPeretto etail 120061 
120071) and NGC 1333 dWalsh et al.|[2006l 120071) This suggests 
that protoclusters often start their evolution from "col d", out-of- 
equilibrium initial conditions (cf. lAdams et al.l 12006]). possibly 
as a r esult of the influence of external triggers (cf. iNutter et ahl 
2006). 



6. Summary and conclusions 

We carried out a detailed observational study of the kine- 
matics of the L1688 protocluster condensations in the central 
Ophiuchus molecular cloud using N2H + (l-0) as a tracer of dense 
gas. Additional observations were also taken in molecular lines 
such as H 13 CO + (1-0), DCO + (2-l), HCO + (3-2), CS(2-1), and 
C 34 S(2-1). Our main results and conclusions are as follows: 

1. The N2H + (l-0) line was positively detected toward 41 of the 
57 compact starless condensations identified by MAN98 at 
1.2 mm, as well as 3 Class 0/Class I protostars (VLA 1623, 
CRBR 85, YLW 15/IRS 43). The same objects were also 
detected in H 13 CO + (1-0) and/or DCO + (2-l) when observed 



in these transitions. In addition, the Class I sources LFAM 26 
and GY210 were detected in H 13 CO + (1-0). 

2. For 29 of the 44 condensations/protostars detected in 
N2H + (l-0) we are confident that the line emission is at least 
partly associated with the compact 1 .2 mm continuum object 
as opposed to the more extended parent cloud/DCO + core, 
since the sources remained positively detected after subtrac- 
tion of the local background line emission. Furthermore, 
at least 17 starless condensations and 3 protostars were 
found to have well-defined line counterparts in position- 
velocity space and could be identified in our background- 
subtracted N 2 H + (101-012), H 13 CO + (1-0), or DCO + (2-l) 
da ta cubes using the Gaussclumps clump-finding algorithm 
of lStutzki & Gustenl ( 1 1990b . 

3. The measured N2H + (l-0) linewidths are narrow, indicative 
of small nonthermal velocity dispersions. The condensations 
of Oph A, Oph Bl, Oph C, Oph E, and Oph F are character- 
ized by subsonic levels of internal turbulence (o~nt/o~t < 1), 
while the condensations of Oph B2 have at most "transonic" 
internal turbulence (cr NT /cr T < 2). 

4. We detected the classical spectroscopic signature of in- 
fall motions toward six 1.2 mm continuum condensations 
(A-SM2, B2-MM16, C-MM5, C-MM6, E-MM2d, and E- 
MM4). For these objects, the optically thick CS(2-1), CS(3- 
2), H 2 CO(2i2 - In), and/or HCO + (3-2) lines are double- 
peaked with a stronger blue peak, while the low-optical- 
depth N2H + (101-012) line peaks in the dip of the optically 
thick line. The same signature was tentatively observed in 10 
other condensations. 

5. The virial masses derived from our N2H + (l-0) detections of 
the condensations generally agree within a factor of ~ 2 with 
the masses estimated from the 1.2 mm dust continuum. On 
this basis, most (~ 75-90%) of the L1688 condensations de- 
tected in N 2 H + (l-0), including essentially all (~ 90%) of 
those more massive than ~ 0.35 M , appear to be gravita- 
tionally bound and prestellar in nature. The status of the star- 
less 1.2 mm condensations less massive than ~ 0.1 M Q is less 
clear since these were often undetected in N2H + . 

6. We estimate that the prestellar condensations of L1688 are 
characterized by a range of lifetimes between ~ 2 x 10 4 yr 
and ~ 5 x 10 s yr. There is however no systematic depen- 
dence of lifetime on mass, so that the observed mass spec- 
tru m is not severely a ffected by the timescale bias pointed out 
bv lClark et al.l (12007b . In particular, the steep, Salpeter-like 
slope of the condensation mass distribution at the high-mass 
end appears to be robust. 

7. Based on the observed distribution of N2H + (l-0) line cen- 
troid velocities, a global velocity dispersion ct\d < 0.4 
km s _1 was estimated for the condensations of L1688. This 
condensation-to-condensation velocity dispersion is sub- 
virial. It implies a crossing time t cross ~ 1 Myr for the con- 
densations within the parent protocluster and a typical colli- 
sion time t co u ~ 1 — 10 Myr between condensations. Since 
these timescales are longer than the estimated condensation 
lifetime ~ 0.2-5 xlO 5 yr, we conclude that, in general, the 
L1688 prestellar condensations do not have time to orbit 
through their parent DCO + core and interact with one an- 
other before evolving into PMS objects. 

8. Using an estimated mean relative velocity of ~ 0.3 - 
0.4 km s _1 between protostellar envelopes and the back- 
ground ambient gas, the mass accretion rate expected from 
competitive accretion was found to be significantly smaller 
than that resulting from gravitational collapse at the Class 
stage. On the other hand, the typical tidal-lobe radius of the 
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L1688 condensations/protostars was estimated to be compa- 
rable to the actual radius of the condensations and protostel- 
lar envelopes as measured in the dust continuum at 1.2 mm. 
This suggests that the overall gravitational potential of the 
protocluster does play an important role in limiting the con- 
densation/envelope masses, but that competitive accretion at 
the protostellar stage is not the dominant mechanism respon- 
sible for determining the final masses of stellar systems in the 
central Ophiuchus star-forming cloud. 
9. We find that competitive, Bondi-like accretion is more likely 
to operate at the prestellar stage and propose that it may 
partly govern the growth of starless, self-gravitating conden- 
sations initially produced by gravoturbulent fragmentation 
toward an IMF, Salpeter-like mass spectrum. 
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